Walks on Apollonian networks 
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We carry out comparative studies of random walks on deterministic Apollonian networks (DANs) and random 
Apollonian networks (RANs). We perform computer simulations for the mean first passage time, the average 
return time, the mean-square displacement, and the network coverage for unrestricted random walk. The diffu- 
sions both on DANs and RANs are proved to be sublinear. The search efficiency for walks with various strategies 
and the influence of the topology of underlying networks on the dynamics of walks are discussed. Contrary to 
one's intuition, it is shown that the self-avoiding random walk, which has been verified as an optimal strategy 
for searching on scale-free and small-world networks, is not the best strategy for the DAN in the thermodynamic 
limit. 
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In the past few years, much attention has been devoted to 
the characterization and modelling of a wide range of complex 
systems that can be described as networks (jj 0|3| • The topo- 
logical properties of real-world networks have been studied 
extensively. But an even more intriguing task, and a natural 
extension of these studies, is to understand how the topologi- 
cal structure of networks affects dynamics taking place on top 
of them 1 4]. Many dynamic processes have been studied on 
complex networks, such as epidemic spreading [5], percola- 
tion 1 6], and synchronization |7] et al. These researches show 
that topologies of networks play a crucial role in determining 
the system dynamical features. 

Random walk has been used for modelling various dy- 
namics in physical, biological, and social contexts |8[. It 
could also be a mechanism of transport and search on net- 
work I9l. ll (Ah ill . Thus, one interesting problem is to study the 
dynamical behavior of a random walker on networks with dif- 
ferent topological properties. Much is known about random 
walks on both regular and random networks 1121 El. Re- 
cently, there have also been several studies of random walks 
on small-world networks (SWNs) fl4l IBl O O and on 
scale-free networks (SFNs) 191 u&L Il9t l20ll . The impacts of 
the heterogeneous topological structures of the networks on 
the nature of the diffusive and relaxation dynamics of the ran- 
dom walk were probed j^Q^ElQ- Also, finding efficient 
algorithms for searching on networks with different topolog- 
ica l prope rties is an important issue related to random walks 

In this paper we investigate walking processes taking place 
upon the determinist Apollonian network (DAN) |23] and its 
variation, the random Apollonian network (RAN) 1 24] . The 
DAN can be defined based of the ancient problem of filling 
space with spheres, first tackled by the Greek mathematician 
Apollonius of Perga [25]. That is, starting with an initial array 
of touching disks, which have curvilinear-triangle interstices, 
disks are added inside each existed interstice in the present 
configuration, such that these disks touch each of the disks 
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bounding the curvilinear triangles. These added disks give 
rise to three smaller interstices, which will be filled in the 
next generation. This process is then repeated for successive 
generations. The DAN is constructed based on this process 
by considering each disk as a node, and the disks in contact 
as the corresponding nodes connected. For each new node 
added to a certain triangle (corresponding to the curvilinear- 
triangle interstice), the three vertices are linked to, and three 
new triangles are created in the network, into which nodes will 
be inserted in the next generation. Different from recursive 
constructing of the DAN 1 23], the RAN starts with a triangle 
containing three nodes. Then, at each time step, one trian- 
gle is randomly selected to add a new node 1 24] . Both net- 
works are simultaneously scale-free, small-world, Euclidean, 
and space filling, and have attracted an increasing interest re- 
cently 0E1 ii. 

We carry out the random walk along the bonds of a given 
network as follows: (i) There is only one walker on the net- 
work at a time; (ii) The random walker is injected onto a 
randomly chosen node on the network, a new node for each 
walker. We will call this node the "origin" of the walk; (iii) At 
each discrete time step t, the walker will jump to one nearest 
neighbor of its current node according to one certain strategy; 
(iv) We average over different random walkers and realiza- 
tions of the network until the results converge. 

The walk strategies adopted by the walker include the fol- 
lowing: random walk (RW), no-back (NB) walk, no-triangle- 
loop (NTL) walk, no-quadrangle-loop (NQL) walk, and self- 
avoiding (SA) random walk. For the RW, the walker may un- 
restrictedly jump to a neighbor node by randomly taking one 
of the links. It forgets all information about its past. The NB 
walk implies that a random walker, if possible, will not return 
to the node it was situated at the previous step. Similarly, NTL 
and NQL random walks mean that the walker will try to avoid 
walking in loops, with three or four edges, respectively, un- 
less there is no other choice. We mention that the NQL walk 
also includes the NTL, which means it eliminates quadrangle 
loops as well as triangle loops. Finally, the SA random walk 
implies that the walker is the smarter. It tries to avoid revisit- 
ing the node that it has already visited in a run of the search. 
Surely, the SA walk also includes the NTL and NQL walks. 
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FIG. 1: (Color online) The mean first passage times (Ty) (upper 
panels) and the average return time (Tu) (lower panels), both of 
which are average over the results of given fc's, for DANs (left-hand) 
and RANs (right-hand). Fitted linear relations of average return time 
(solid lines) are obtained with a slope of —1 for both DANs and 
RANs. 

One can imagine that an unrestricted RW walker may be 
trapped into high-clustering regions and revisit nodes there 
for a long time by reason of walking in loops. While a walker 
adopting other strategies like NB, NTL and NQL can escape 
from those regions easier, and therefore improve the search 
efficiency. Thus, different clustering properties of the under- 
lying networks will induce different efficiencies of the above 
mentioned strategies. In return, the clustering effect of the 
networks could be reflected by the efficiencies of those search 
strategies. 

In the context of transport and search, the mean first passage 
time (MFPT) and the average return time (ART) are important 
characteristics of the random walk. The average return time 
denoted by (Tu) is the average time needed for a walker to 
return to the origin i, which is a special case of MFPT from 
i to j denoted by (Ty). Noh and Rieger 1191 have proved 
that the MFPT is negatively correlated with Kj, and (Tu) is 
determined only by the total number of links and the degree 
of the node i 

(Tn)=M/Ki (1) 

with A/" = U = 1,...,N). Namely, nodes with 

higher degrees are visited earlier and more frequently, and 
therefore targets on these nodes are more easily found out than 
on nodes with smaller degrees. In Fig. ^ we depict simulations 
of the MFPT (the upper panels) and ART (the lower panels) 
as a function of nodes degree for the DAN and RAN with size 
N = 9844. After averaging over the MFPT and ART for a 
given degree k, the linear property of the log-log plot of ART 
(the lower panels) is in excellent agreement with Eq. Q, and 
the negative correlation between MFPT and k are also pre- 
sented (the upper panels). 

In the following, we will investigate walking processes 
on DANs and RANs. The corresponding results on Watts- 
Strogatz (WS) (K = 3, p = 0.1) [29] and Barabasi-Albert 



TABLE I: The PSP threshold P c and average path length (L) of the 
four networks with size TV = 9844. 
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(BA) (to = 77io = 3) J30ll networks are also presented for 
comparison. Note that, the four networks have the same aver- 
age degree (k) = 6. 

Firstly, let us discuss the search efficiency of all the five 
walk strategies mentioned above for the four networks. We 
suppose at every step, the walker adopting various strategies 
only know neighbors of its present node. So if the target is 
at one neighbor of the present node where the walker stays, 
this round of search is over. The walks for the five strategies 
are performed and the search times (see Fig. [2} are averaged 
over 200 different couples of initial positions and targets of 
the walker, and 50 realizations of the networks. One can find 
some common properties from Fig. [2] For example, the ap- 
proximate power-law relations hold for all search processes 
of RW, NB, NTL, NQL and SA walks. As another example, 
the RW walk is the most inefficient search strategy for all the 
four networks. Nevertheless, some obvious differences still 
exist among them. The search efficiencies of NB, NTL and 
NQL walks nearly collapse into one only for BA network (see 
Fig. 0(c), also in Ref. |20[). This can be understood by its 
comparatively smallest clustering coefficient as to the other 
three networks. In the weak clustering networks, the NB, 
NTL and NQL walkers, who adopt increasingly strict law to 
avoid revisit, can not take their advantages and their search ef- 
ficiencies keep very closely or even collapse into one. Hence, 
for DANs and RANs, the efficiency-improvement of NTL and 
NQL walks as to the NB walk reflects the fact that high clus- 
tered nodes are popular on them. 

The SA walk was proved to be generally the most efficient 
strategy for searching on networks, whose small-world prop- 
erties were considered to be the reason 1 20] . However, it is 
counterexample for the DAN (see Fig. |2jb)). The DAN is one 
of typical regular networks which has deterministic size when 
the number of the generations is certain 11231 12611 . Its short 
average length and high clustering coefficient indeed describe 
a small-world scenario. Fig. |2jb) shows the average search 
times for DANs with 3 to 8 generations, corresponding to net- 
work sizes 16, 43, 124, 367, 1096, and 3283, respectively. 
Contrast to the conclusions in Ref. 1 20], the SA walk does not 
reduce the search time and is even more inefficient than the 
NTL and NQL walks when is large. We argue that what is 
important for the improvement of the search efficiency is not 
taking more strict law to eliminate repeating visits, but con- 
sidering the factual structure of the networks. 

According to the former analysis of MFPT, we know that 
the high-degree nodes would be visited early, so as to be 
avoided by SA walker preferentially. Thus, to some extent the 
SA search process can be considered as the question of inten- 
tional attack 13111 or preferential site percolation (PSP) on net- 
works. The intentional attack means the removal of nodes and 
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FIG. 2: (Color online) Average search times of various search strategies for RANs (a), DANs (b), BA networks (c), and WS networks (d). 
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FIG. 3: Mean-square distance (R 2 ) as a function of time for DANs, 
RANs, WS and BA networks with size TV = 9844 and average de- 
gree (k) — 6 fixed. The dot line with slope 1 is plotted for compari- 
son. 



their incident edges aiming at those greatest degree nodes, and 
the PSP threshold P c is a measure of how robust the network is 
against this attack. When the network size is relatively small 
(N < 10 5 ), DANs and RANs are much frailer than SW and 
SF networks under the intentional attacks (see Table. |l]for P c ), 




FIG. 4: (Color online) Network coverage N co v after t steps on (solid 
line from top to bottom) BA networks, DANs, RANs and WS net- 
works with N = 9844 and (k) = 6 fixed. 



because of the crucial importance of the high-degree nodes to 
the network integrates 1 24] . When walking on networks, if the 
walkers strictly avoid visiting these crucial high-degree nodes, 
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they can hardly move from their current part of the network 
to other parts, as if the network was intentionally attacked 
and disintegrated into partes. However, as a search strategy, 
the SA search processes on the original network are not to- 
tally the same as the intentional attacks, because SA walker 
at node i can revisit nodes, including these crucial nodes, al- 
though, with probability 1/ki only in condition that all the ki 
neighbors have already been visited. The compulsive revisits 
insure that the targets in different parts of networks can finally 
be found. Considering the similarity between SA search and 
intentional attack or PSP, we can conclude that the SA search 
strategy, compared with other mentioned strategies, will im- 
prove the search efficiency to certain extent for the networks 
with higher PSP threshold. 

Next, we will present our simulation results for two quan- 
tities, the mean-square displacement and the network cover- 
age (the average number of distinct sites visited). The mean- 
square displacement (R 2 (t)) of a particle diffusing in a given 
space, which is a measure of the distance R covered by a typ- 
ical RW walker after performing t steps, is one of the most 
basic quantities in the random walk theory lllc! fl9ll . In most 
cases, this quantity is described by an expression of the form 
(R 2 (t)) ~ t a . The value of the parameter a classifies the 
type of diffusion into normal linear diffusion (a = 1), sub- 
diffusion (a > 1), or superlinear diffusion (a < 1). When 
considering the unrestricted RW walk along the bonds of net- 
works, the maximum value of a can be 2 I19B . To calculate 
(R 2 (t)) we first, at each time step, find the minimal distance 
from the current position of the walker to the origin (i.e., the 
smallest number of steps needed for the walker to reach the 
origin) using a breadth-first search method. Then we allow 
the walker to move through the network until (R 2 (t)) has sat- 
urated. Finally, the results are average over different initial 
positions of the walkers and realizations of the network. We 
simulate (R 2 (t)) for the DAN, RAN, WS and BA networks 
with N = 9844, and report it as a function of time in Fig. 
[5] One important feature is the fact that {R 2 (t)) equilibrates 
after a few steps to a constant displacement value. This is 
a simple manifestation of the small diameter of these finite 



networks. Note also that, because of the differences of the 
average path length (L) of these networks (see Table. [J for 
(£)), the plateau values are also different. For the DANs and 
RANs, one can find that the slope of (R 2 ) in the early time is 
small, especially that of RANs. This can be explained as a re- 
sult of their highly clustering effect which induces the walker 
initially spends many time exploring the cluster it was created 
in and thus the distance to the origin increases slowly. Af- 
ter a few steps the walker escapes its initial territory and dif- 
fuses around the entire network. From the slope of (R 2 ), we 
know that diffusions both on DANs and RANs are sublinear, 
and diffusions both on BA and WS networks are superlinear. 
The network coverage denoted by N cov is defined as average 
number of distinct visited nodes of the RW walker. Results 
of N cov on the four networks are presented in Fig. |4] The 
number of steps performed is nearly two orders of magnitude 
smaller than the size of the networks, in order to decrease fi- 
nite size effects. The clustering effects of DANs and RANs 
are reflected by their relatively low coverage compared to BA 
networks. 

In summary, we present comparative studies of the dynam- 
ics of random walks on deterministic Apollonian networks 
and random Apollonian networks with other networks. The 
efficiency of search strategies are simulated and the clustering 
and short path effects of networks to the dynamics are dis- 
cussed. It is found that, for the DAN with large size, the best 
search strategy is no longer the self-avoiding random walk. 
The importance of the high-degree nodes to the network inte- 
grates are suppose to be one reason. Thus the search strategy 
on networks should be assigned closely based on the topolog- 
ical structure of the underlying networks, including the clus- 
tering, the average path length as well as the essentiality of 
the high-degree nodes. Since search is a problem of extreme 
importance for so many natural and artificial networks, this 
finding may be valuable in practicability. Finally, the simu- 
lation results of mean-square displacement (R 2 (t)) and net- 
work coverage N cov also show the influences of the structure 
of networks to the dynamics of random walks. 
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